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Abstract: 

To match the precision of present and future measurements of Z-boson produc- 
tion at hadron colhders, electroweak radiative corrections must be included in the 
theory predictions. In this paper we consider their effect on the transverse mo- 
mentum (pt) distribution of Z bosons, with emphasis on large px- We evaluate, 
analytically and numerically, the full one-loop corrections for the parton scatter- 
ing reaction qq Zg and its crossed variants. In addition we derive compact 
approximate expressions which are valid in the high-energy region, where the weak 
corrections are strongly enhanced by logarithms of s/M^. These expressions include 
quadratic and single logarithms as well as those terms that are not logarithmically 
enhanced. This approximation, which confirms and extends earlier results obtained 
to next-to-leading logarithmic accuracy, permits to reproduce the exact one-loop 
corrections with high precision. Numerical results are presented for proton-proton 
and proton-antiproton collisions. The corrections are negative and their size in- 
creases with pt- For the Tevatron they amount up to —7% at 300 GeV. For the 
LHC, where transverse momenta of 2 TeV or more can be reached, corrections up 
to —40% are observed. We also include the dominant two-loop effects of up to 8% 
in our final LHC predictions. 
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1 Introduction 



The study of gauge boson production has been among the primary goals of hadron 
coUiders, starting with the discovery of the W and Z bosons more than two decades 
ago [l]. The investigation of the production dynamics, strictly predicted by the 
electroweak theory, constitutes one of the important tests of the Standard Model. 
Furthermore, this reaction contributes to the background for many signals of new 
physics. Being embedded in the environment of hadronic collisions, the reaction 
necessarily involves hadronic physics, like parton distributions, and depends on the 
strong coupling constant. In turn, the cross section for the production of W and Z 
bosons and their rapidity distribution can be used to gauge the parton distribution 
functions [2 which are important ingredients for the prediction of numerous reaction 
rates. 

Differential distributions of gauge bosons, in rapidity as well as in transverse mo- 
mentum (pt), have always been the subject of theoretical and experimental studies. 
In the region of small px, multiple gluon emission plays an important role and contri- 
butions of arbitrary many gluons must be resummed to arrive at a reliable prediction 
OH]. At larger transverse momenta the final state of the leading order process con- 
sists of a or Z boson plus one recoiling jet. QCD corrections in next-to-leading 
order to this process are mandatory for the correct description and can amount to 
several tens of per cent depending on the observable under consideration, including 
jet definition, as well as the renormalization and factorization scales [Sj. The eval- 
uation of next-to-next-to-leading order corrections involves two-loop virtual plus a 
variety of combined virtual plus real corrections and is a topic presently pursued by 
various groups (see e.g. Ref. [H]). 

For the experiments at the Large Hadron Collider (LHC) a new aspect comes into 
play. The high center-of-mass energy in combination with the enormous luminosity 
will allow to explore parton-parton scattering up to energies of several TeV and 
correspondingly production of gauge bosons with transverse momenta up to 2 TeV 
or even beyond. In this region electroweak corrections from virtual weak boson ex- 
change increase strongly, with the dominant terms in L-loop approximation being 
leading logarithms of the form log^^(,s/M^), next-to-leading logarithms of the 
form a^\og^^~\s/M^), and so on. These corrections, also known as electroweak 
Sudakov logarithms, may well amount to several tens of percent. They have been 
studied in great detail for processes involving fermions in Refs. [H |H]. Investigations 
on the dominant and the next-to-leading logarithmic terms are also available for 
reactions involving gauge and Higgs bosons [HIEJIEI- A recent survey of the liter- 
ature on logarithmic electroweak corrections can be found in Ref. The impact 
of these corrections at hadron colliders has been studied in Ref. (T^I. Specifically, 
hadronic Z-boson production at large has been investigated in next-to-leading 
logarithmic approximation, including the two- loop terms Numerical results for 
the complete one-loop terms have been presented in Ref. [T5] . 

It is the aim of this work to obtain an independent evaluation of the complete 
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one-loop weak corrections to the same reaction, and to present the full result in 
analytic form. At the partonic level the reactions qq ^ Zg, qg ^ Zq and qg — > Zq 
with q = u,d,s,c or b have to be considered which are, however, trivially related by 
crossing and appropriate exchange of coupling constants. We split the corrections 
into an "Abelian" and a "non-Abelian" component. The ultraviolet divergences of 
the former are removed by the renormalization of the fermion wave functions as 
expected in Abelian theories. The latter receives also divergent contributions from 
the renormalization of the Z-boson wave function and the electroweak parameters. 
The electroweak coupling constants are renormalized in the framework of the mod- 
ified minimal subtraction (MS) scheme. We present analytic results for the exact 
one-loop corrections that permit to predict separately the various quark-helicity 
contributions. We also derive compact analytic expressions for the high-energy be- 
haviour of the corrections. Here we include quadratic and linear logarithms as well 
as those terms that are not logarithmically enhanced at high energies but neglect 
all contributions of 0{M^/s). The accuracy of this approximation is discussed in 
detail. 

After convolution with parton distribution functions, radiatively corrected pre- 
dictions for transverse momentum distributions of Z bosons at hadron colliders are 
obtained. Concerning perturbative QCD, these predictions are based on the lowest 
order and thus proportional to as- To obtain realistic cross sections, higher-order 
QCD corrections would have to be included, in next-to- leading or even next-to- next- 
to-leading order. 

The paper is organized as follows: In Sect. |21the Born approximation, our con- 
ventions, the kinematics, and the parton distributions are introduced. Sect. 01 is 
concerned with a detailed description of the radiative corrections. In Sect. 13. ll the 
strategy of our calculation is described and the relevant Feynman diagrams are intro- 
duced. Sect. 13. 21 is concerned with the algebraic reduction of the amplitudes to a set 
of basic Lorentz and Dirac structures that are multiplied by scalar form factors and 
scalar one-loop integrals. Subsequently the renormalization is performed (Sect. I3.3|l 
and a compact form for the results is given (Sect. 1^3), with a decomposition into 
Abelian and non-Abelian contributions. Special attention is paid to the behaviour in 
the high-energy limit (Sect. 13. 5|) . which can be cast into a compact form. Indeed the 
quadratic and linear logarithms confirm the evaluation of the Sudakov logarithms 
obtained earlier in Ref. fl]. Sect. EJthen contains a detailed discussion of numerical 
results, both for pp and pp collisions at 14 TeV and 2 TeV, respectively. The lead- 
ing and next-to-leading two-loop logarithmic terms are included in this numerical 
analysis. We also discuss the sensitivity of the results towards the choice of the 
renormalization scheme and justify the use of the MS scheme adopted in this paper. 
Sect. El concludes with a brief summary. 
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2 Conventions and kinematics 

The pt distribution of Z bosons in the reaction h^h^ ^ 2' + jet is given by 

—. = '^^1/ dx2 6'(xiX2 - 'rmm)/fti,i(a;i,/i^)//i2,j(a;2,/i^)-5 — , (1) 

dpT 77 "'0 Jo apT 



where fmin = (pt + '^t)^/s, mx = + is the colhder energy. The 

indices i,j denote initial state partons {q,q,g) and fhi,i{x, ^'^), fh2,j{x, fi"^) are the 
corresponding parton distribution functions, o"'-' is the partonic cross section for the 
subprocess ij Zk and the sum runs over all i,j combinations corresponding to 
the subprocesses 

qq Zg, qq Zg, gq Zq, qg Zq, qg Zq, gq Zq. (2) 

The Mandelstam variables for the subprocess ij — ^ Zk are defined in the standard 
way 

s = {pi+ Pjf, i = {pi - pzf, u = {pj - pzf. (3) 

The momenta pi, pj, pk of the partons are assumed to be massless, whereas p| = M|. 
In terms of Xi,X2,Pt and the collider energy ^/s we have 

s = XiX2S, t = — (1-COS6'), u=— (1 + COS6'), (4) 



with cos 6* = Y 1 — 4p^s/(s — M^Y corresponding to the cosine of the angle between 
the momenta pi and pz in the partonic center-of-mass frame. 

The angular and the distribution for the unpolarized partonic subprocess 
ij — >■ Zk read 

^^i^^M"? (5) 



dcos0 327rA/:ys2 



and 



da*-' Pt 



dpT 9,'KNijS\t - u 



(6) 



where 



E = ^EE (7) 

pol col 

involves the sum over polarization and color as well as the average factor 1/4 for 
initial-state polarization. The factor l/A^jj in (0)-®, with Nqq = Ngg = N^, Ngq = 
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Nqg = Nqg = Ngq = Nc{N^ — 1), and Nc = 3, accounts for the initial-state colour 
average. The factor 



i-u = (s- M|) cos 6* = s^J {xiX2 - rmin)(a;iX2 - r^m + 4pt"^t/s) (8) 

in the denominator of (jH)) gives rise to the so-called Jacobian peak, which arises at 
cos 6* = or, equivalently, at xiX2 = Tmin and is smeared by the integration over the 
parton distribution functions. 

The unpolarized squared matrix elements for the processes (j21) are related by the 
crossing-symmetry relations 



9<?|2 ^ _ 



qq\2 



99|2 ^ _ 



El-^ 



qq\2 



and 



j2\M'f = El-^^'f 

Moreover, as a result of CP symmetry. 

Using these symmetries we can write 
da^i^2 



t<r-*U 



(9) 
(10) 



dpi 



X 



+ 



V [ dxi [ dx2 9{xiX2 - f„ 
^„„Jo Jo 



q=u,d,c,s,b ' 

.2\ r / ..2\ 



)fh2,q{x2,^^ ) + (1 ^ 2) 



d(T^"5 

dpT 



fhuq{Xl,fi^)fh2,g{^2, /i^) + fhuqi^l, /^^ ) A2 ,9 (3^2 , /i^) +(1^2) 



By means of and is expressed in terms of X^I-^^'^P) i-e. the explicit 

computation of the unpolarized squared matrix element needs to be performed only 
for one of the six processes in 0. 

To lowest order in a and as, for the qq Zg process 



A=L,R 



tu 



(13) 



where a = e^/ (47r) and as = 9s/ (4^) are the electromagnetic and the strong coupling 
constants and represents the eigenvalue of the generator associated with the Z 
boson in the representation corresponding to right-handed (A = R) or left-handed 
(A = L) quarks q\. In terms of the weak isospin T^^ and the weak hypercharge Yq^ 




with the shorthands Cw = cos^w and Sw = sin6'w for the weak mixing angle 

In general, for gauge couplings we adopt the conventions of Ref . [12] . With this 
notation the gqq vertex and the Vqq vertices with V = A,Z, read 

-i^s^"7^ Vwv = ieY E ^x^L',^ (15) 

g /" Q' / A=R,L 

where ux are the chiral projectors 

'^R = ^(l+75), t^L = ^(l-75), (16) 

are the Gell-Mann matrices and are matrices in the weak isospin space. For 
diagonal matrices such as we write /^^/^ = 5qqil^^. The triple gauge-bosons 
vertices read 




3; 



+ g^^^^^^{ks-k,r% (17) 



The definition of the antisymmetric tensor e as well as useful group-theoretical 
identities can be found in App. B of Ref. [TU] . 



3 0{a) corrections 

In this section we present the one-loop weak corrections to the process qq Zg. 
The algebraic reduction to standard matrix elements and scalar integrals is described 
in Sect. 13.21 Renormalization, in the MS and the on-shell scheme, is discussed in 
Sect. 13.31 In Sect. 13.41 we present analytical results for the unpolarized squared 
matrix elements, and the high-energy behaviour of the corrections is derived in 
Sect. HSl 



3.1 Preliminaries 

As discussed in the previous section, the 6 different processes relevant for Z + 1 jet 
production are related by crossing symmetries. It is thus sufficient to consider only 
one of these processes. In the following we derive the one-loop corrections for the 
qq — > Zg process. The matrix element 

Mf = Ml'' + 6Mf (18) 

is expressed as a function of the Mandelstam invariants 

S={Pq+Pqf, i={Pq-pzf, U= {Pq-Pzf- (19) 
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■fWWW 



-isiasLajisiy 



(tl) 



vwvw^ 



(t2) 



Figure 1: Tree-level Feynman diagrams for the process qq Zg. 



The Born contribution, A^q'^, results from the t- and u-channel diagrams of Fig. d 
The loop and counterterm (CT) diagrams contributing to the corrections, 

SMf = 6MZ^^^^ + 6Ml%^, (20) 

are depicted in Fig. |21 and Fig. El respectively. 

We do not include electromagnetic corrections, i.e. we restrict ourselves to the 
virtual weak contributions of (9(a^as). In the kinematical region where the Z-boson 
transverse momentum is non-vanishing the emission of the gluon does not give rise 
to soft or coUinear singularities and the (renormalized) virtual weak corrections are 
finite. All (real and virtual) quarks appearing in the diagrams of Fig. El are treated 
as massless^ and diagrams involving couplings of quarks to Higgs bosons or would- 
be-Goldstone bosons are not considered. Quark-mixing effects are also neglected. 
The only quark-mass effects that we take into account are the mt-and mfe-terms in 
the gauge-boson self-energies, which contribute to the counterterms. 

Our calculation has been performed at the level of matrix elements and provides 
full control over polarization effects. However, at this level, the analytical expressions 
are too large to be published. Explicit results will thus be presented only for the 
unpolarized squared matrix elements 

El-^fr = El-^o'r + 2Re [E(-^o"T 6MT]+Oia'as). (21) 



3.2 Algebraic reduction 

The matrix element ()18j) has the general form 

Mf = iegst^ v{p,)M\'^''u^u{p,)el{pz)e*M), (22) 

A=R,L 

where the 7^-terms are isolated in the chiral projectors uj\ defined in Since 

""^ Quark- mass contributions are of 0(TO^/Af|) and can thus be neglected for q ^ t. The con- 
tributions of 0{m^/M'^) arising from those diagrams of Fig. |21 that involve initial-state h quarks 
and virtual W bosons are large at the partonic level. However, at the hadronic level, also these 
contributions can be neglected since the subprocesses initiated by b quarks are suppressed by the 
very small &-quark density of protons. 
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Figure 2: One-loop Feynman diagrams for the process qq —>■ Zg. The diagrams 
v5, v6 and b3 involve only charged weak bosons, W^, whereas the other diagrams 
receive contributions from neutral and charged weak bosons, V = Z, W^. 




(cl) 




(c2) 



AAAAAA/ 



c3) 



www^ 



(c4) 



<XlAAAAA/ 



->JLQ_Q_SJ2J2^ 

c5) 



-(K^AAA/V^ 

(c6) 



Figure 3: Counterterm diagrams for the process qq Zg. 



we treat quarks as massless, M.^''^'^ consists of terms involving an odd number of 
matrices 7^* with p = 0, . . . , 3. The polarization dependence of the quark spinors 
and gauge-boson polarization vectors is implicitly understood. 
In analogy to (fTHj) . (plj) we write 

Ml''"' = M^'''" + SM^''"' = ^A^t'ioops + '^-^tcT- (23) 

The unpolarized squared matrix element ()2H1 is obtained from TVIq '^'^ and (5A1^''^^ 
using 

J {MI'Y 5Mr = n-'aasiNl - 1) E Tr [^.TWo'^V.-^-^i'"''^'" 

A=R,L 

x«„.L.-»V (24) 



Ml 



The Born contribution reads 



-'^'0 — -'^a'-'O ' '-'0 — f \ 7 ) {'^'^) 

and is simply proportional to the gauge group generator J^, which determines the 
dependence of Aio on the weak mixing angle, the chirality (A = R, L) and the 
weak isospin of quarks {u,d). Let us now consider the combinations of gauge group 
generators that appear in the loop diagrams of Fig. |21 The diagrams involving 
virtual Z bosons are simply proportional to (I^)^- Instead, the diagrams with virtual 
W bosons (diagrams sl-b3 with V = W^) yield combinations of non-commuting 
gauge-group generators and triple gauge-boson couplings W^z^ which 

can be expressed as^ 

jzjw^jw^ ^ jw^jw^jz (diagrams si, s2, vl, v2), (26) 
^^w^w^zjw^jw^ ^ (diagrams v5, v6, b3), (27) 

jw^jzjw^ ^ jZjW^jWT _ c^^3 (diagrams v3, v4, bl, b2), (28) 

where fermionic indices of the generators as well as summations over are 
implicitly understood. The contribution of the loop diagrams of Fig. |21can thus be 
expressed as a linear combination of and T^, 



6 Ml 



loops - ^ 



v=z,w± 



(29) 



^The relation (|2ZI) can be derived from Eqs. (B.5) and (B.24) in Rcf. 10 whereas (jSHl follows 
from (|27(l combined with the commutation relation [I^ ,1^] = — (i/s^)e^ w^ZjW _ 
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with 

The amphtudes 6A'^\^^{My) are denoted as the Abehan (A) and non-Abehan (N) 
contributions, since the latter originates from the non-commutativity of weak inter- 
actions, whereas the former is present also in Abelian theories. 

These amplitudes have been reduced algebraically using the Dirac equation, 
the identity p^s^{p) = for gauge-boson polarization vectors and Dirac algebra. 
Moreover, tensor loop integrals have been reduced to scalar ones by means of the 
Passarino-Veltman technique. The result has been expressed in the form 

10 14 

s-^'^aMm^) = T.T.n/M)sruMi). (31) 

i=l jr=0 

where the form factors J-^^i^{My) are rational functions of Mandelstam invariants 
and masses. The tensors 





= l^{Tj>z - 


-#.-)7^ 




= 






= YPz, 




5r 


= -7X^ 






= 1%, 






= -7X. 






= i:^z-i 






= i.iz-1 






= {iz-i 




'->10 


= {iz-i 


K)v%fz. 



correspond to the massless subset of the standard matrix elements of Ref. ^B] and, 
apart from J^iMy) = 1, the functions Jj{My) represent the scalar one- loop integrals 
resulting from the reduction.^ The one- and two-point functions, 

Ji(M^) = Ao(M2), 
J2(M2) =i?o(M|;0,0), 
J3(M2) =So(M|;M2,M2), 
MM^) =5o(s;0,0), 
J5(M2) =5o(m;M2,0), 

MM^) = Bo{i;M^,0), (33) 

■^For the scalar integrals Aq, Bq, Cq and Dq we adopt the notation of FeynCalc ^21- However, 
we choose their normalization according to Ref. [16) . i.e. we include the factor (27r/i)'*~^ which is 
omitted in the conventions of FeynCalc. 
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are ultraviolet (UV) divergent. The resulting poles, in D 
amount to 



where Sq^ is the Born amplitude defined in ((251) • These singularities are cancelled 
by the counterterms (see Sect. 13. 3|) . The remaining loop integrals are free from UV 
singularities. The three-point functions 



uv 



:'->o ) 



2e dimensions, 
(34) 



Co(S,0,0;0,0,M2), 

Co(m,M|,0;M2,0,0), 

Co(n,M|,0;0,M2,M2 

Co(t,M|,0;M2,0,0), 

Co(t,M|,0;0,M2,M2) 



(35) 



are finite for non-vanishing Z-boson transverse momentum. In addition, we obtain 
the three- and four-point functions 



Co(m,0,0;M2,0,0), 
Co(t,0,0;M2,0,0), 
Co(S,M|,0; 0,0,0), 



Do(0,0,M|,0,m,s;M2, 0,0,0), 
Do(0,0,M|,0,t,S;M2, 0,0,0), 
Do(M|,0,0,0,t,M;M2,M2,0,0), 



(36) 



which result from the reduction of the box diagrams bl-b3. As a consequence of 
vanishing quark masses, these latter integrals ()3(i|l give rise to mass singularities of 
collinear nature. Such singularities originate from the propagators of those massless 
(virtual) quarks that couple to the massless (real) gluon. 



(37) 




specifically from the integration region with xp'^, where the momenta of these 
quarks become collinear to the gluon momentum. However, the box diagrams are 
finite since the quark-gluon vertex (jH7|l yields 



XiX — 



gfgifg 



x{x - 1) 2^g{p, 



0, 



(38) 



in the collinear limit. Indeed, as a consequence of cancellations between the singu- 
larities from the Cq and the Dq functions in ()36|) . our result is finite. 

In order to control these cancellations at the analytical and numerical level, we 
have regulated the singularities by means of an infinitesimal quark-mass parameter 
A. Then, using Ref. [IHI we have expressed the singular parts of -Do functions through 
singular Cq functions and checked that these cancel against the singular Cq functions 
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resulting from the reduction. Finally we have written our result in terms of finite 
combinations of Dq and Co functions, 



2 n ^ ..2 n ^2 n^ " M|)Co(w, M|, 0; , A^, 0) 



SU+ [t + U)My 

uCo{u, 0, 0; M^, A^ 0) + (g - M|)Co(g, M|, 0; 0, 0, A^) 
su + (t + m)M^ 

, . 9 9 ,9 ^ (t-M2)Co(^',M2 0;M2,A2,0) 
Ji3(M2 ) = Do(0, 0, M|, 0, t, s- Ml, 0, A^, 0) - ^ /; , ^ ^2 

tCo(t, 0, 0; M^, A^ 0) + (g - M|)Co(g, M|, 0; 0, 0, A^) 
si+{t + u)Ml 
MM^) = Do(M|, 0, 0, 0, t, n; , , 0, A^) 

iCo{i, 0, 0; M^, 0, A^) + «Co(m, 0, 0; , A^, 0) 



iu - {i+u)M, 



(39) 



smg 



which correspond to the finite parts of four-point integrals expressed as Dq — Dq 
with the singular parts D^^^ in terms of Cq functions. It was checked that the 
functions J12, J13, J14 are indeed finite and numerically stable for X/Mz -C 1. The 
numerical results presented in Sect. IU have been obtained using X/Mz = 10~^. 



3.3 Renormalization 

The renormalization of the process qq —>■ Zg is provided by the diagrams depicted 
in Fig. ini The counterterms (CTs) that are responsible for the contributions of 
diagrams cl, c2, c3 and c4 read 



i^ J2 ^x^Z, 

A=R,L 




E 

A=R,L 



(40) 

Since there is no (9(a)-contribution to the renormalization of the strong coupling 
constant gs, these CTs depend only on the fermionic wave- function renormalization 
constants SZg^ [see fl43|) ]. Their combined contribution to the qq Zg process, 
i.e. the sum of diagrams cl, c2, c3 and c4, vanishes. The renormalization of the 
qq — > Zg process is thus provided by the diagrams c5, c6, which originate from the 
Zqq counterterm. 




and yield 



ier E 

A=R,L 



Qx 



6Z. 



Qx 



(41) 



+ 



^Qx^'-'Qx 



'^0 5 



(42) 
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with iSq"^ defined as in (j2SI)- As indicated in (jH]), the Zqq CT can be expressed as a 
hnear combination of and T^^, analogous to ()29j) . The corresponding coefficients 
involve the wave-function renormalization constants of massless chiral quarks,^ 



6Z, 



-Re 



S«'^(0) 



a ^ 

y=z,w± 



9a 



Ml 



(43) 



where only weak corrections {V = Z, W^) are included, and the following combina- 
tions of renormalization constants 



6C, 



SC. 



9a 



The renormalization constants associated with the Z-boson wave function read 



(44) 



5Z 



zz 



-Re 



5Z 



AZ 



-2Re 



(M|) 



(45) 



and have been evaluated using the self-energies of Ref . [inj . 

For the renormalization of e and Cw we have considered two different schemes, 
the MS and the on-shell scheme. In the MS scheme, the CTs for e and read 



5c^ 


MS 


a 






47r 




MS 




e2 




4:71 



22s2 

^ + 2(p 



uv 




where 



l/e - 7e + log(47r) + log 



(46) 



(47) 



in D = 4 — 2£ dimensions and p = M^/(c^M|). In (P7|) we have included a loga- 
rithmic term that renders the renormalized amplitude independent of the scale /ir 
of dimensional regularization. This is equivalent to the choice /zr = Mz within 
the usual MS scheme. We also note that the counterterms ()46|1 contain contribu- 
tions proportional to {p — 1) which depend on the relation between M^/ and Mz- 
In principle, in an 0{a) MS calculation the value of the W mass, which appears 

'^'^'^'■^{p'^) are the chiral components of the one-particle irreducible two-point function of massless 
fermions, 



1+ J2 ^aS'^^p') 



A=R,L 
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only in loop diagrams, should be derived from the on-shell Z-boson mass and the 
MS mixing angle using the tree-level relation Mw = c^^Mz- In this case the (p — 1) 
contributions in PU)) would vanish. However, it seems more natural to use the 
on-shell value of My/ as an input parameter in our calculation. This violates the 
tree-level mass relation and introduces loop corrections that are implicitly contained 
in the numerical value of Mw- Their effect on our 0{a) predictions is formally of 
C(a^) since Mw does not contribute at tree level. This procedure is thus consistent 
at 0{q). However, it affects also those 1/e poles that appear in 5Zaz since they 
depend on Mw- Such singularities give rise to terms proportional to (p — l)/e that 
are again, formally, of 0{a^) and are compensated by the (p — 1) terms that we 
have introduced in P6|) . This procedure removes all higher-order effects from the 
divergent part of the counterterms (jSj), which becomes independent of (p — 1), and 
ensures the cancellation of the singularities (jH^ originating from the loop diagrams 
for arbitrary values of p. 

In the on-shell (OS) scheme, the weak mixing angle is defined as = M^/M| 
and the corresponding CT reads 



5^ 



w OS 



Re 



T.^^{M^ 



(4J 



As on-shell input parameter for the electromagnetic coupling constant we used a 
a{Mz), defined as 



a(M: 



a(0) 



1 - Aa(M|)' 



Aa(M: 



Re 



(49) 



where a(0) = e(0)^/(47r) is the fine-structure constant in the Thompson limit and 
nfem represents the fermionic contribution to the photonic vacuum polarization. 
The CT associated with = 47ra(M|) is given by 

2sw S4^(0) 



5e^ OS MO)' ^^,,,2. 



e{Of 



Re 



nt^(Mi) +ii^^M- 



M\ 



2 ' 



(50) 



where 5e(0) is the CT in the Thompson limit, and Hj^^(O) represents the bosonic 
contribution to the photon propagator. 

Independently of the renormalization scheme, the above CTs yield the ultraviolet 
singularities^ 



5Z, 



a 1 

A-K e 



y=Z,W± 



5C, 



UV 



0, sc, 



a 



UV 



2713% € 



(51) 



that cancel the 1/e poles ()34j) originating from the loop diagrams. The scheme 
dependence of the complete result is discussed in Sect. HI 

^Note that SC^^ is free from ultraviolet singularities, i.e. the pole associated with the Abelian 
coupling structure in 1)42(1 results exclusively from the fermionic wave-function renormalization 
constant 6Zq^ whereas the poles originating from the renormalization of the electroweak coupling 
constants {Se, (5cw) and the Z-boson field renormalization constants {SZzz, SZaz) cancel in 6C^^. 



13 



3.4 Result 



Let us present the complete 0{a'^as) result for the unpolarized squared matrix 
element fj21|) for the qq — > Zg process. Combining the contributions ()29|). 
and using we obtain 



a 1 



X 



A=R,L 



v=z,w± 



+ 



27r si 



w) 



where the coupling factors are specified by and (jSHI), and 

Pzt^PzA _ t' + u' + 2sMl 



Ho = Itt l^qSo" :^gS^'''' 



9v 



iu 



(52) 



(53) 



represents the Born contribution. The counterterms dC^^^ are defined in ()44|) - 

dSOD. i/f^^(M^) are the Abelian (A) and non-Abelian (N) contributions resulting 
from the loop diagrams of Fig. |21 and the fermionic wave-function renormalization 
constants (EHl). These functions can be written as 



Ht^^M^) = Re E [</^(M^) + <?:,(M^)] J,(M^, 



(54) 



where Jo (My) = 1 and the remaining Jj{My) are the loop integrals defined in (j33|) . 
fl35j) . (jnni)- The fermionic wave- function renormalization constants yield [cf. 



M2 



(55) 



and K^Yii^v) — otherwise. The coefficients resulting from the loop integrals of 
Fig. El read 



1 10 

oEnV(M^)Tr 



9/111' 



i=l 



PZ/iPZ/i' 

Ml 



(56) 



where jF^yp^(My) are the form factors of (jSH). The explicit expressions for these 
coefficients are presented in Appendix El 

In fact, as a consequence of Ward identities (see Sect. I3.(i)l . the contributions of 
the pz/iPz/i'/M^ terms in (jKHjl and are equal zero. We note that the contribu- 
tions of right- and left-handed quarks to the unpolarized squared matrix element 
can be easily recognized as the terms with A = R and A = L, respectively. 
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3.5 High-energy limit 

In this section we discuss the behaviour of the corrections to the qq — > Zg 
process in the high-px region. More specifically, we derive approximate expressions 
for the functions H^^^ and the counterterms SC^J^ for the case where all energy 
scales are much higher than the weak-boson mass scale, i.e. in the limit where 
My^/s — > and t/s, u/s are constant. 

In this approximation, which is applicable for transverse momenta of order of 
100 GeV and beyond, the one-loop weak corrections are strongly enhanced by loga- 
rithms of the form log(s/M^). The logarithmically enhanced part of the corrections, 
which constitutes the so-called next-to-leading logarithmic (NLL) approximation, 
was already presented in Ref. This NLL part consists of double- and single- 
logarithmic terms that result from the functions H^^^ in (jK^ and read^ |14j 



log' 



log' 




(57) 



whereas the counterterms do not contribute to the NLL approximation, i.e. 



r^A NLL e^N NLL p. 



(5J 



The NLL corrections ()57|1 are proportional to the Born contribution Hq defined in 
(j^Hjl . The above result is free from logarithms associated with the running of the 
electroweak coupling constants. Also no single logarithms originating from virtual 
collinear weak bosons that couple to the final-state Z boson appear in ()57|) . This is 
a general feature of gauge-boson production processes [01 HD] • 

In the following, we provide the complete asymptotic behaviour of the weak 
corrections, which includes, in addition to double and single logarithms, also those 
contributions that are not logarithmically enhanced in the high-energy limit. This 
constitutes formally the next-to-next-to-leading logarithmic (NNLL) approximation 
and provides a precision of order (M^/s) log'(s/M^). 

Let us first consider the functions Hi^^ in (|52|) . which result from the loop 
diagrams of Fig. |21and the fermionic wave-function renormalization. The asymptotic 
behaviour of these functions was derived using the general results of Ref. jTH] for 
the high-energy limit of one-loop integrals. In addition, in order to simplify the 



^The correspondence between (|5^ . (|^ and (5), (10) of Ref. ^l] can be easily seen by means 
of the relation 

y=z,w± 



9A 



where C™ 



is the electroweak Casimir operator and the —Q^ term represents the subtraction of 
the electromagnetic contributions Ji]. We also recall that the result of Ref. is based on the 
equal-mass approximation, Mz ~ Myy- 
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dependence of the result on the masses of the Z and W bosons, we have performed 
an expansion in the Z- W mass difference using = 1 — / M| as the expansion 
parameter and keeping only terms up to the first order in s%j. 

The resulting NNLL expressions for iff (M^), E^{Ml/) and B^^M^) have the 
general form 



H^I''{Ml) ^1^^ Re 



A/N 



{Ml\ 



til 



u 



tu 



+ 9^"" {Ml) 



(59) 



i.e. they involve the function + v?)/iu^ which corresponds to the high-energy 
limit of the Born contribution Hq [see fl53|) ]. and two other rational functions of the 
Mandelstam invariants, which describe different angular dependences. The functions 
gf^^ consist of logarithms of the kinematical variables and constants. For the non- 
Abelian contribution H^{M^r) we obtain 



g^{M^) = 2Auv + lo; 



W) 



w) 




log' 
'73 



-u 



+ 4, 



log" ( ^1 +log" ( Ij +log(T] +log(| 



log' 



All' 



(60) 



where Auv is defined in (H7jl . and H^^M^) is absent in (fH^ . For the Abelian 
contributions H^{My) with V = Z,W we obtain 



9^iMl 



gtiM'y] 



log' 



+ 3 log 



+ log^j +log[| 



-9^{M'w) 



+ 



log 




(61) 



The results (j60|) - ()6Hl . for the qq Zg process, are valid for arbitrary values of the 
Mandelstam invariants. They can thus be easily translated to all other processes 
in (j21) by means of appropriate permutations of s, t and m, which are specified by 
the crossing relations (jH)) and ()1U|). We note that the arguments of the logarithms 
in (jUm) -()6H ) are in general not positive. Logarithms with negative arguments are 
defined through the usual ie prescription, f — »• f + ie for f = s,t, u. In practice we 
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need only the real parts of the logarithms, 



Re log (a;) 
Re log^(x) 



log(|x|), 

log'(|x|) -7r2^(-x). 



(62) 



The TT^-terms originating from the analytic continuation of the double logarithms 
depend on the signs of their arguments and these latter are determined by the signs of 
the Mandelstam invariants. For the qq — > Zg process, one has s > 0, t < and m < 0. 
However, these signs change when the Mandelstam invariants are permuted in order 
to translate the results (jUnj) - (jM|) for the qq —>■ Zg process to the other partonic 
processes in The vr^ terms resulting from (jU^ are thus process dependent. 

Simple explicit expression that apply to all processes in (0) can be obtained using 
the relations 



is) + e{i) + e{u) 



9{s)9{i) = e{i)e{u) = e{u)e{s) 



0, 



(63) 



which are valid for s > 0, t < and -u < as well as for any permutation of the 
Mandelstam invariants (crossing transformations). Combining ()6()|l - (pTH|) we obtain 



Re^o''(^^) = 2Auv + log' 



Re^f(M^ 
Reg^iM, 



w) 



log' 



log' 



\u\ 



+ log' 



\u\ 



3 
2 

log' 



log' 



+ log' 



\u\ 



\u\ 



-log' 



log' 



+ 



TX 



-|4 + 9»M)]-^ + 4. 



e{i) - e{u) 



(64) 



and 



Reg^{Ml 



RegtiMl) 
Reg^iM^r) 



-log' 



" V31og^ " 



Ml + 2 



+ logf|^j +log(H 



log' 



+ log' 



^ [2 + 3^^(-s)] - I 



-Ref7f(M^) + ^ 



-Reg^{M, 



w) 



\u\ 
\s\ 



(65) 



For the qq Zg and qq Zg processes, with s > 0, t < and m < 0, all 6 
functions in ()64|l - (jHK|l vanish. The ^-terms become relevant when the above results 
are translated to the processes that involve gluons in the initial state. 
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Let us compare the double and single logarithms of My that appear in (jMjl - 
dnni) with the NIL result of Ref. [Hj, which is given in In the M^/s 

limit the two results agree up to contributions resulting from the difference between 
(yf|^(M^) and Qq^M"^). Since the calculation of Ref. was based on the equal-mass 
approximation, My/ = Mz, we conclude that the two results are consistent. 

Let us now consider the contribution of the counterterms (CTs) 6C^J^ in (j52|) . 
These CTs [see consist of on-shell gauge-boson self-energies and are in- 

dependent of the high-energy scales s, t and u. They are thus free from large loga- 
rithms. Since the exact analytical expressions for the CTs are relatively complicated 
as compared to the compact NNLL contributions we present simple ap- 

proximate expressions for 6C^J^ . Restricting ourselves to the case of the MS scheme, 
we use: (i) A heavy-top expansion including terms of order \og{mt) and (l/mt)°; 
(ii) A first-order expansion in the Z-W mass difference, using 1 — M^r/M"^ = as 
expansion parameter; (iii) The light-Higgs approximation Mh = Mz- The result- 
ing expressions for the renormalization constants associated with the Z-boson wave 
function read 



6Z 



AZ 



a 



47r S\mC\> 



7 + 34s,2 



+ 44(p - 1) A 



+ (17 - 36s^ 



'Ts ~ 27 



(78 - 89s 



6Zzz = ^ 



a 1 



5-104 + 46st^ 



uv 



84) log 



64 



uv 



6 



log 



+ 



11-374 ^ 



V3 



113 - 5734 
36 



(66) 



And for the CTs 



within the MS scheme we obtain 



MS a 1 



^^A MS « 



15 



477 Sl, 

113 



+ 



36 

2534 TT 



49s^ 

-log 



mi 



105 4567 



12 



8 



+ 



216 



5C. 



N MS 




-Auv + 



2o2 p2 



84) log 



:i7-36s^ 



TT 



^(78 -894) 



mi 



(67) 



As one can easily infer from the numerical values^ SZ^z = —2.316 x 10 ^, SZzz = 
3.943 X 10-3, 6Zaz = -2.888 x lO'^, SZzz = 4.533 x lO'^, the contribution of the 
CTs to the cross section is of order of a few permille whereas the error associated 
with the approximations ()66|) . (jF)7j) is below the one permille level. 

These numerical values are obtained using the MS input parameters given in Sect. 01 and 
Auv = 0. 
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3.6 Checks 



In order to control the correctness of our results we performed various consistency 
checks. We have verified that the one-loop corrections ()31|) satisfy the Ward Identi- 
ties 



PZi^PguViPg) [5^5',A/n(^V-)^a] U{pg) 



0, 



0. 



(6J 



Similar Ward identities hold for the lowest-order amplitude. The cancellation of the 
ultraviolet divergences and the compensation of the fictitious coUinear singularities 
discussed in Sect. 13.21 have been verified analytically and numerically. The high- 
energy limit of our result has been worked out analytically and we have shown that 
the leading- and next-to-leading logarithmic contributions agree with Ref. |14j . 

Moreover, the calculation has been performed in two completely independent 
ways, based on different computer-algebra and numerical tools. On one hand we 
have written algorithms for the algebraic reduction in Mathematica [201 and used a 
set of routines provided by A. Denner for the numerical evaluation of loop integrals. 
On the other hand we have performed the reduction by means of FeynCalc [12j 
and used FF [21 to evaluate loop integrals. The results have been implemented 
in different Fortran codes and comparing them we find very good agreement at 
numerical level. 



4 Predictions for the hadronic Z production at 
high transverse momentum 

In the following we discuss numerical predictions for the quantities calculated above. 
The lowest order (LO) and the next-to-leading-order (NLO) predictions result from 
(HSI) and respectively. The next-to-leading- logarithmic (NIL) 

approximation, at one loop, corresponds^ to (j37j) - (j^ . The next-to-next-to- leading- 
logarithmic (NNLL) predictions, also at one loop, are based on (j^ - (jU7j) . All results 
are obtained using LO MRST2001 parton distribution functions (PDFs) j22j. We 
choose p^ as the factorization scale and, similarly as the scale at which the running 
strong coupling constant is evaluated. We also adopt the value as(^I) = 0.13 and 
use the one-loop running expression for asifJ^"^), in accordance with the LO PDF 
extraction method of the MRST collaboration. We use the following values of pa- 
rameters Mz = 91.19 GeV, Mw = 80.39 GeV, = 176.9 GeV, nib = 4.3 GeV, 
Mh = 120 GeV. For the calculation in the MS scheme we use a = 1/128.1, 

®For details concerning the treatment of angular-dependent logarithms at the NLL level we 
refer to Ref. 
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Figure 4: Relative one-loop corrections to the partonic differential cross sections 
do"*-' /d cos 6* at cos 6* = for (a) uu channel, (b) dd channel, (c) gu channel, (d) gd 
channel. The sohd, dashed and dot-dashed lines denote the modulus of the TZ ratios, 
as defined in the text, for the full NLO cross section, the NLL approximation and 
the NNLL approximation of the one-loop cross section, respectively. 



= 0.2314, whereas for the on-shell scheme a = 1/128.9, = 1 — = 
1 — M^^/M"^ are taken. The NLO predictions use the form of the counterterms 
given by ()44j) - (jHn|) . whereas the NNLL results are based on the approximate ex- 
pressions (inni)-(jnZ|) for the counterterms. Unless otherwise noted, the results are 
obtained using the MS scheme. 

We begin by investigating NLO, NLL and NNLL relative corrections to the 
partonic (unpolarized) differential cross section do"*-'/dcos6' [see ©]. To this end 
we define 

dq-NLo/dcosg 
d(3"Lo/dcos6' 

and similarly T^nll/nlo "^nnll/nlo- These ratios, calculated at cos6' = 0, 
are displayed as a function of \^ in Fig. El We consider four processes: uu^Zg 
(Fig. it), dd^Zg (Fig.ib), gu~^Zu (Fig.ih), gd^Zd (Fig. Hi). The size of the 
full weak NLO correction grows with the energy and reaches 30% for all channels 
at \fl = 4 TeV. The sign of this correction is negative. From Fig. |3 we conclude 
that the NLL terms provide a fairly good approximation to the full NLO result for 
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> 200 GeV, with the remaining terms responsible for less than 3% of the cross 
section in the uu and dd channels, and less than 1% in the gu and gd channels. 
The quality of the NNLL approximation is very good in all channels, better or 
comparable to 1% in the full region under consideration. The absolute LO and 
NLO cross section and the TZ^^ ratios for specific values of cos^, and different 
collision channels are listed in Table H 

The transverse momentum distribution da/dp^ at the LHC is shown in Fig. 
We display separately the absolute values of the LO, NLO, NLL and NNLL differ- 
ential cross sections (Fig. E^) and the relative correction wrt. the LO result for the 
NLO, NLL and NNLL distributions (Fig. Eb). The relative correction wrt. the LO 
is now defined as 

-TP had _ dq-NLo/dpT -I 

- daLo/dpT " ^ ^ 

for the NLO case, and similarly for the NLL and the NNLL cross sections. The 
quality of the NLL and NNLL high-energy approximations is shown in more de- 
tail in Fig. Et. The importance of the NLO corrections increases significantly with 
Pt- The NLO correction results in a negative contribution ranging from —13% at 
Pt = 500 GeV up to —37% at = 2 TeV of the LO cross section. We observe that 
the NLL approximation works very well. It differs from the full NLO prediction by 
about 1% at low p^ and by 0.2% at ~ 2 TeV, cf. Fig Et- The quality of the 
NNLL approximation is at the permille level (or better) in the entire p^ range. 

In view of the large one-loop effects, we include the dominant two-loop terms 
|14j . In Fig. iniwe show the relative size of the corrections in the NLO approximation 
(^NLO/LO' solid line), and in the approximation which includes the next-to- leading 
logarithmic two- loop terms ("T^nnlo/lo' dotted line). These additional two- loop 
terms are positive and amount to 8% for = 2 TeV. 

To underline the relevance of these effects, in Fig. [7| we present the relative 
NLO and NNLO corrections for the cross section, integrated over p^ starting from 
Pt = Pt"*, as a function of p^p*. This is compared with the statistical error, defined 
as AcTstat/c" = with N = Cx BR(Z— >/, i/^) x ctlo- We include all leptonic 

decay modes of Z, corresponding to BR = 30.6%, and assume a total integrated 
luminosity C = 300 fb~^ for the LHC [21]. It is clear from Fig. [71 that the size of the 
one- and two-loop corrections is much bigger than and comparable to the statistical 
error, respectively. 

We also perform a similar analysis for high transverse momentum Z production 
at the Tevatron. In Fig. |Hlwe show the transverse momentum distribution (Fig. |H^), 
the relative size of the corrections (Fig. |H|d) and the quality of the one-loop NLL 
and NNLL approximations (Fig. |Ht). In Fig. IHlthe relative size of the corrections to 
the integrated cross section is compared with the statistical error expected for an 
integrated luminosity C = llfb~^ At the energies of the Fermilab collider, the 
NLO weak correction is of the order of the statistical error and should be taken into 
account when considering precision measurements. The two-loop terms turn out to 
be negligible. 
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da'^/dcose [pb] 


X 103 




GcV 


cos 9 


LO 


NLO 


NLO/LO 


NLL/NLO 


NNLL/NLO 


uu 


500 


0.0 


1.7012 X 10° 


1.6341 X 10° 


-39.431 


19.501 


9.7186 




1000 


0.0 


3.5934 X IQ-i 


3.2261 X 10-1 


-102.22 


17.397 


6.6197 




2000 


0.0 


8.1081 X 10-2 


6.5185 X 10-2 


-196.05 


17.518 


4.1422 


uu 


500 


0.5 


2.8203 X 10° 


2.7210 X 10° 


-35.224 


25.906 


7.2284 




1000 


0.5 


6.0587 X 10-1 


5.4862 X 10-1 


-94.497 


24.664 


4.9420 




2000 


0.5 


1.3718 X 10-1 


1.1185 X 10-1 


-184.66 


25.881 


3.2372 


dd 


500 


0.0 


2.1930 X 10° 


2.1031 X 10° 


-40.978 


19.759 


11.038 




1000 


0.0 


4.6322 X 10-1 


4.1349 X 10-1 


-107.36 


16.723 


7.4723 




2000 


0.0 


1.0452 X 10-1 


8.2852 X 10-2 


-207.31 


16.002 


4.5613 


dd 


500 


0.5 


3.6357 X 10° 


3.5013 X 10° 


-36.970 


25.829 


8.3255 




1000 


0.5 


7.8103 X 10-1 


7.0278 X 10-1 


-100.18 


23.838 


5.6448 




2000 


0.5 


1.7684 X 10-1 


1.4205 X 10-1 


-196.74 


24.376 


3.5721 


gu 


500 


0.0 


6.9404 X 10-1 


6.5844 X 10-1 


-51.284 


-6.9286 


-1.0724 




1000 


0.0 


1.6267 X 10-1 


1.4171 X 10-1 


-128.83 


-8.2907 


-0.2447 




2000 


0.0 


3.7676 X 10-2 


2.8750 X 10-2 


-236.92 


-9.8681 


0.7254 


gu 


500 


0.5 


5.9307 X 10-1 


5.5148 X 10-1 


-70.127 


-13.649 


0.9367 




1000 


0.5 


1.3822 X 10-1 


1.1654 X 10-1 


-156.82 


-16.198 


1.0469 




2000 


0.5 


3.1936 X 10-2 


2.3182 X 10-2 


-274.09 


-19.650 


1.4449 


gd 


500 


0.0 


8.9468 X 10-1 


8.4822 X 10-1 


-51.931 


-0.0567 


-2.0723 




1000 


0.0 


2.0969 X 10-1 


1.8206 X 10-1 


-131.76 


-0.6210 


-0.7120 




2000 


0.0 


1.8569 X 10-- 


3.6717 X 10-- 


-211.01 


-1.3163 


0.1704 


gd 


500 


0.5 


7.6453 X 10-1 


7.1421 X 10-1 


-65.819 


-5.1329 


-0.1256 




1000 


0.5 


1.7818 X 10-1 


1.5103 X 10-1 


-152.38 


-6.6209 


0.5373 




2000 


0.5 


4.1168 X 10-2 


2.9994 X 10-2 


-271.41 


-8.7165 


1.1609 



Table 1: Absolute value of the LO and NLO partonic differential cross section 
d(T*-'/dcos6' and the ratios "T^nlo/lo' "^nll/nlo ^'^'^ "^nnll/nlo in permille, for the 
partonic processes uu ^ Zg, dd ^ Zg, gu — > Zu and gd Zd. The running strong 
coupling q;s(a*2) is taken at the scale /j,'^ = = {1 — cos^ 6){1 — Mz/sfs/A. 
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Figure 5: Transverse momentum distribution for pp^Zj at ^/s — 14TeV. (a) 

LO (solid), NLO (dashed), NLL (dotted) and NNLL (dot-dashed) predictions, (b) 
Relative NLO (solid), NLL (dotted) and NNLL (dot-dashed) weak correction wrt. 
the LO distribution, (c) NLL (dotted) and NNLL (dot-dashed) approximations 
compared to the full NLO result. 
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Figure 6: Relative NLO (solid) and NNLO (dotted) corrections to the px distribu- 
tion for pp^Zj at ^/s = 14 TeV. 




Figure 7: Relative NLO (solid) and NNLO (dotted) corrections wrt. the LO predic- 
tion and statistical error (shaded area) for the unpolarized integrated cross section 
for pp^Zj at = 14 TeV as a function of p^* . 
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Figure 8: Transverse momentum distribution for pp^Zj at \/i — 2 TeV. (a) LO 
(solid), NLO (dashed), NLL (dotted) and NNLL (dot-dashed) predictions, (b) Rela- 
tive NLO (solid), NLL (dotted) and NNLL (dot-dashed) weak correction wrt. the LO 
distribution, (c) NLL (dotted) and NNLL (dot-dashed) approximations compared 
to the full NLO result. 
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Figure 9: Relative NLO (solid) and NNLO (dotted) corrections wrt. the LO and 
statistical error (shaded area) for the unpolarized integrated cross section for pp^Zj 
at = 2 TeV as a function of p^*. 

Finally, in Fig. El we illustrate the dependence of the NLO px distribution on 
the choice of the renormalization scheme. At low the results in the on-shell and 
MS schemes differ by around 1% and the difference grows with reaching 6% at 

= 2 TeV. This effect is mainly due to the different treatment of the weak mixing 
angle in the two renormalization schemes. As well known, the relation between the 
MS and on-shell definitions of the weak mixing angle is provided by the p parameter 

as 

^ = mk = P^ S^-l = ^4 = 5^Ap, (71) 

where p = (1 — Ap)~^ and the symbols with and without hat denote MS and on-shell 
quantities, respectively. The input parameters used in our calculation, = 0.2314 
and = 1 — M^/M| ~ 0.2228, correspond to As^ ~ 3.8% and are extracted 
from precision electroweak measurements taking all available loop corrections into 
account. Instead, loop corrections beyond 0{a) are not included in our calculation 
and, in particular, the deviation observed in Fig. is due to missing two-loop 
(and higher-order) corrections related to the p parameter. The scheme dependence 
resulting from a/s"^ terms amounts to 

where A'^^^s^, = {c^/ s'^)Ap^^^ ~ 4.5% corresponds to the one-loop corrections to the 
p parameter, which are included in our on-shell predictions through the counterterm 
(HHJ. This scheme dependence (f7^ is thus due to the higher-order contributions 
As^ — A^^^^s^ ~ —0.7%. Their relatively large size results from the combined effect 
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Figure 10: Relative difference of the NLO transverse momentum distribution for 
pp^Zj calculated in the MS and on-shell (OS) schemes at i/i = 14 TeV. 

of 0{aasm'f) j2lj, 0{aalmf) jSH], 0{a'^mf) ^ and 0{a^m'^t) jHO] corrections to 
the p parameter and is consistent with the effect observed in Fig. El at small px- 
In addition, the scheme-dependence resulting from the one-loop logarithmic terms 
is of order 




This effect is due to missing two-loop corrections of order Ap alog^(s/M^). Its size 
is proportional to As^ ~ 3.8% and grows with energy. This explains the high-px 
behaviour in Fig. ^| 

We stress that the effects (f7^ - ()73p are entirely due to missing higher-order 
terms related to Ap and that such missing terms concern only the calculation in 
the on-shell scheme. Indeed, Ap enters our predictions only through the relation 
between the weak mixing angle and the weak-boson masses in the on-shell scheme 
whereas the MS calculation does not receive any contribution from Ap. The large 
scheme-dependence in Fig. ^Jhas thus to be interpreted as large uncertainty of the 
one-loop prediction in the on-shell scheme whereas such uncertainties are absent in 
the MS scheme. This motivates the choice of the MS scheme adopted in this paper. 

5 Summary 

In this work we have calculated the one-loop weak corrections to hadronic production 
of Z bosons at large transverse momenta. Analytical results are presented for the 
parton subprocess qq Zg and its crossed versions. Special attention has been de- 
voted to the high-energy region, s 3> M^, where the weak corrections are enhanced 
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by logarithms of s/M^r. For this region we have derived approximate expressions 
that include all large logarithms as well as those terms that are not logarithmically 
enhanced. The quadratic and linear logarithms confirm earlier results obtained to 
next-to-leading logarithmic accuracy. The complete high-energy approximation pre- 
sented in this paper is in excellent agreement with the exact one-loop predictions. 
We give numerical results for proton-antiproton collisions at 2 TeV (Tevatron) and 
proton-proton collisions at 14 TeV (LHC) in the region of large tranverse momentum 
(Pt)- The corrections are negative and their size increases with px- At the Tevatron, 
transverse momenta up to 300 GeV will be explored and the weak corrections may 
reach up to —7%. At the LHC, transverse momenta of 2 TeV or more are within 
the reach. In this region the corrections are large, —30% up to —40%, and even the 
dominant two-loop logarithmic terms must be included in any realistic prediction. 



A Analytical results 



Here we present the explicit analytical expressions for the coefficients Kf^^ {My), 
which are defined in ()56p and correspond to the Abelian (A) and non-Abelian (N) 
contributions resulting from the Feynman diagrams of Fig. |21 
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Non-Abelian coefficients 
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